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Abstract 

Vibronic effects in resonant electron transport through single-molecule junctions are analyzed. 
The study is based on generic models for molecular junctions, which include electronic states on the 
molecular bridge that are vibrationally coupled and exhibit Coulomb interaction. The transport 
calculations employ a master equation approach. The results, obtained for a series of models with 
increasing complexity, show a multitude of interesting transport phenomena, including vibrational 
excitation, rectification, negative differential resistance (NDR) as well as local cooling. While some 
of these phenomena have been observed or proposed before, the present analysis extends previous 
studies and allows a more detailed understanding of the underlying transport mechanisms. In 
particular, it is shown that many of the observed phenomena can only be explained if electron- 
hole pair creation processes at the molecule-lead interface are taken into account. Furthermore, 
vibronic effects in sytems with multiple electronic states and their role for the stability of molecular 
junctions are analyzed. 

PACS numbers: 73.23.-b,85.65.+h,71.38.-k 



I. INTRODUCTION 



For more than a decade, molecular electronics [IHS] has been a very active and chal- 
lenging field of research. One of the basic ideas is to exploit the diversity of molecules and 
the possibilities of modern synthesis to design molecular systems with specific functions 
for nanoscale electronic devices. Another motivation is the possibility to investigate single 
molecules under controllable nonequilibrium conditions. Various techniques, including me- 
chanically controlled break junctions [T0l - [T5] . electro- migrated molecular junctions [T6T42T] . 
scanning tunneling microscopy [22H27], and very recently, on- wire lithography in combina- 
tion with in-situ 'click chemistry' [28, 29J, have been employed to contact a single molecule 
with two macroscopic electrodes. Once such a molecular junction is established, external 
electric fields, either a bias or a gate voltage [15], (TTJ ETJ ED] , can be used to investigate the 
conductance of a single molecule. A variety of interesting transport phenomena have been 
proposed and experimentally observed [31 |5] EH9l EI], including, e.g., switching behavior 
[SMS], rectification [TQl EHH ES] and negative differential resistance [JJJ HH1 EB EZHSH] • 

However, a detailed understanding of the experimental results, especially in the reso- 
nant transport regime, where electrons may populate an intermediate state on the bridging 
molecule, has not been achieved yet. One of the interesting and challenging aspects of 
electron transport in molecular junctions is the intricate interplay between the electronic 
and vibrational degrees of freedom. Because of the small size of molecules, the charging of 
the molecular bridge is often accompanied by significant changes of the nuclear geometry 
that indicate strong coupling between electronic and vibrational degrees of freedom. As a 
consequence, the vibrational modes of a molecular junction can be highly excited result- 
ing in significant nonequilibrium effects [lOTHS"] . This aspect distinguishes single-molecule 
junctions from traditional quantum dot systems [HH17]. Vibrational signatures indicating 
strong vibronic coupling as well as strong excitation of vibrational modes were identified for 
a number of molecular junctions [JBJ HBrtZU 1221 123 ESI Novel experimental tech- 

niques based on measuring the force needed to break a junction [M] or employing Raman 
spectroscopy [551 EE] allow a characterization of the current-induced vibrational nonequilib- 
rium state of a single-molecule junction. These data complement the information carried by 
the respective current- volt age characteristics. 

Various theoretical approaches have been employed to describe vibrationally coupled 
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electron transport through single molecules. While scattering theory approaches [57H60] 
can be used to address the regime of strong molecule-lead coupling, nonequilibrium Green's 
function approaches [KH HTl loTHBT] additionally allow a non-perturbative description of the 
associated nonequilibrium state of such a junction, especially with respect to the vibrational 
degrees of freedom. Numerically exact methods, based on path integrals [6E1 EH] or multi- 
configurational wave-function methods [70], provide valuable insights and benchmarks for 
specific model systems and problems that may not be addressed by perturbation theory or 
other approximative schemes. Master equation approaches [H2J [7THE2] , although perturba- 
tive with respect to the coupling between the molecule and the leads, have been proven to 
be very powerful, as they are capable of describing all interactions on the molecular bridge 
accurately and, for simple model systems, very efficiently In this work, we employ a master 
equation approach that is based on a second order expansion in the molecule-lead coupling 
[62] 172] [81] . Master equation approaches that take into account higher-order effects with 
respect to this coupling have already been employed [TU EHHH01 E2]- The corresponding 
higher-order transport processes [101 183H58] , however, are beyond the scope of this work, 
where we focus on the resonant transport regime. 

In this article, we analyze nonequilibrium transport phenomena induced by electronic- 
vibrational coupling in molecular junctions. To this end, we consider generic models for 
vibrationally coupled electron transport including electronic states on the molecular bridge 
that are vibrationally coupled and exhibit Coulomb interaction. The results show a multi- 
tude of interesting phenomena that extend previous studies and may facilitate the interpre- 
tation of experimental results. In particular, it is found that the current-induced vibrational 
excitation in transport through a molecular bridge with a single electronic state increases 
significantly with increasing bias voltage and/or decreasing electronic- vibrational coupling. 
We show that this phenomenon is caused by electron-hole pair creation processes [89H9T] . 
which to the best of our knowledge have not been considered in detail in this context before. 
Further analysis shows that electron-hole pair creation processes can also explain the strong 
enhancement of vibrational rectification effects in situations where the vibrational degree 
of freedom acquires a highly excited nonequilibrium state. This complements the analy- 
sis of the phenomenon of vibrational rectification given in Ref. [6TJ Electronic-vibrational 
coupling may also cause vibrationally induced NDR effects [731 EEJ E21 EH]- In this work, 
we reinterpret the NDR-mechanism outlined in Refs. [TJH in terms of pair creation pro- 
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cesses. Furthermore, a novel mechanism for vibrationally induced NDR is discussed, which, 
in contrast to earlier studies EH ES] , extends over a broad range of bias voltages. 

In molecular junctions, where multiple electronic states participate in the charge trans- 
port, a number of additional vibronic processes have to be considered. In particular, as we 
have shown recently [41 j . higher- lying electronic states facilitate resonant absorption pro- 
cesses that may deexcite the vibrational degrees of freedom. This mechanism reduces the 
current-induced vibrational excitation of a molecular junction and results in local cooling 
[HJ H2]. In the present paper, we give a detailed analysis of effects due to multiple elec- 
tronic states. In this context, we also show that repulsive Coulomb interactions may further 



enhance the stability of a molecular junction (cf. IIIB3). We refer to this phenomenon as 
'Coulomb Cooling'. Since polyatomic molecules typically include numerous active vibra- 
tional modes and often exhibit multiple closely lying electronic states, these phenomena are 
expected to be of relevance for most molecular junctions. 



The article is organized as follows. In Section II A we introduce the model Hamiltonian 
used to describe electron transport through a single-molecule junction. The derivation of 
the master equation approach and expressions for the observables of interest, in particular, 
current- volt age characteristics and the average vibrational excitation of a molecular junction, 
are outlined in Sees. |II B| and |II C[ Explicit formulas adapted to the specific model systems 
are detailed in appendices [X] and |Bj The role of coherences is elucidated in appendix [Cj 



Sec. Ill comprises numerical results and a discussion of the different transport phenomena. 



Thereby, we consider transport phenomena that involve a single electronic state (Sec. Ill A) 



and two electronic states (Sec. Ill B ) . Besides the basic transport mechanisms that are 



discussed in Sec. HI A 1 and HI B 1 we study vibrationally induced rectification in Sees. 



Ill A 2 and III B 2 as well as vibrationally induced negative differential resistance in Sees. 



Ill A 3 and IIIB4 Vibrational excitation processes in junctions with two electronic states 



and their influence on the stability of a molecular junction are analyzed in Sees. Ill B 1 to 



III B 3 Throughout the article we use units where h = 1. 
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II. THEORY 



A. Model Hamiltonian 

We consider electron transport through a single molecule that is bound to two metal 
leads. Such a molecular junction is described by a set of discrete electronic states, which are 
localized on the molecular bridge (M) and interact with a continuum of electronic states in 
the left (L) and the right (R) lead, respectively. The corresponding model Hamiltonian is 
given by 

H c i = Y e i c l c i + Y e k c l c k (!) 

+ Y u iM°i - 6 i)( c h - 6 j) 

+ Y (Vkiclci + h.c). 

fcSL,R;ieM 

Thereby, e& denote the energies of the lead states with corresponding creation and anni- 
hilation operators c\ and Likewise, is the energy of the ith electronic state on the 
molecular bridge, which is addressed by creation and annihilation operators c] and q. The 
coupling matrix elements characterize the strength of the interaction between the elec- 
tronic states of the molecular bridge and the leads and determine the so-called level-width 
functions V K>ij (e) = 2ix J2keK V k\ v kj^(^ ~ e fc) (K=L,R). 

Additional charging energies, due to Coulomb interactions, are accounted for by Hubbard- 
like electron-electron interaction terms, t4,(cjcj — <5j)(cjcj — 8j). Thereby, the parameters Si 
distinguish states that are occupied (d>» = 1) or unoccupied (5, = 0) in the molecular system 
at equilibrium. While in the present study we use a generic model, the first-principles 
determination of the parameters requires the introduction of a reference system [96J . In the 
present paper, we consider the neutral molecule in its electronic ground state as a reference 
system. As a result, for an electronic state above the Fermi-level (8i = 0), the associated 
energy denotes the energy required to add an electron to the ith electronic state of the 
reference system. For an electronic state below the Fermi- level (Si = 1), e« denotes the 
energy required to remove an electron from state i, in accordance with Koopmans' theorem. 
The Fermi energy of the leads is set to €f = eV. 

Upon transmission of electrons, the molecular bridge may be vibrationally excited. We 
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describe the vibrational degrees of freedom of a molecular junction within the harmonic 
approximation, 

#vib = ^2^l a a j a a a + ^ ^iaQa(c\ci- Si), (2) 

a i£M;a 

where the operator a£, denotes the creation operator of the ath oscillator with frequency 
Q a . The coupling between the electronic and the vibrational degrees of freedom is assumed 
to be linear in both the vibrational displacements Q a = a a + and the electron (or hole) 
densities (cjcj — Si) [HHl EEJ [97] . The respective coupling strengths are denoted by X ia . 
Because we employ the normal modes of the ground-state of the neutral molecule, there is 
no coupling between the electronic states and the normal modes of the molecular junction 
in this state. This is imposed in the electronic-vibrational coupling term by the parameters 
Si. The Hamilton operator of the overall system is given by the sum 

H = H c \ + H vih . (3) 

In the limit of vanishing molecule-lead coupling, — > 0, the Hamiltonian H can be 
diagonalized by the small polaron transformation [101 E21 EH EH] 

H = e s He~ s = H S + H B + #sb, (4) 

h s = ^2 z i c l°i + 5^ fi « a « a " 

i a 

+ U ij(4 c i - 5i)(c)cj - Sj), 
i<j 

k 

H S b = yiiVjejXjclci + h.c), 

ki 

with 



S = ^2-^(clci-Si)(al-a a ), (5) 

ia 

Xi = exp[J^-^(a a -4)]. (6) 

a 

Thereby, we have partitioned the transformed Hamiltonian H in three parts, H = H$ + 



+ Hsr, with H$ representing the molecular bridge, the leads and f/sB describing the 
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interactions between the molecular bridge and the leads. Due to the small-polar on trans- 
formation, there is no explicit electronic- vibrational coupling in H^. However, electronic- 
vibrational coupling appears in the transformed Hamiltonian H at three different places: 

• in the polaron-shifted energies: ^ = q + (25, — 1) ^2 a {)$ a /£l a ), 

• in additional electron-electron interactions, which shift the original electron-electron 
interaction terms: = — 2 ^2 a (Xi a Xj a /Q a ), 

• and in the molecule-lead coupling term if sb that is renormalized by the shift operators 



B. Master equation approach 

Density matrices have been proven to be a powerful tool in describing quantum- 
mechanical systems [99til03j . Once the density matrix g of a given system is known, all 
observables O of that system can be obtained from the trace 

(O) = Tr{gO} = H qO \a) = e»bO ba . (7) 

a ab 

Thereby, the elements of the density matrix are given as Q ab = (a\g\b), where |a) and \b) 
are elements of a complete set of orthonormal basis functions that span the Hilbert space 
of the overall system. Since Eq. ^ is invariant under the small polaron transformation, we 
consider in the following H as the Hamiltonian of the system. 

The time-evolution of a density matrix g(t) is determined by the Liouville - von Neumann 
equation: 

^ = -i\H,g(t)]=-iCg(t), (8) 
g(t) = e- iCt g(0), (9) 

where the initial state at time t = is encoded in the respective density matrix g(0). Here, 
L denotes the Liouville operator C g = [H, g\ . 

To describe an open quantum system such as a single molecule coupled to a reservoir of 
electrons (in the left and the right electrode), it is expedient to employ the reduced density 
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matrix of this system p, which is obtained by taking the trace over the degrees of freedom 
of the reservoirs (or baths (B)), 



p(t) = Tr B {f?(t)}. (10) 
Formally, this can be achieved by the projection operator 

PQ(t) = PBTT B Mt)}=p B p(t), (11) 

and its orthogonal complement Q — 1 — P. Here, p B denotes the thermal equilibrium density 
matrix of the reservoir 

p B = Z'V^, Z = Tr B {e-^ B }. (12) 

Assuming a factorized initial condition g(0) = p(0)pB, the equation of motion for the 
reduced density matrix is given by the Nakajima-Zwanzig equation [1041 1105j 



d_ 
Ft 



Pg(t) = -iPCPg(t)- [ drPjCe- iQCT Q£Pg(t-T), 

Jo 

(13) 



which represents a formally exact equation. For practical applications, approximations are 
required to solve this equation of motion. In the present context, we assume that the 
molecule-lead coupling, Hsb, is weak. Employing a second order expansion in Hsb and the 
condition tr B {-f^sBPB} = 0, the Nakajima-Zwanzig equation can be simplified to the master 
equation [H2J ED E21 HHS] 

|p(t) = -i[H s ,p(t)} (14) 

POO 

- drtr B {[i? S B, p S B(r),p(t)p B ]]}, 
Jo 

with 

H sb (t) = e - l ^ s+ " B)T H SB e i(Jls+llB)T . (15) 



To obtain the time-local master equation described by Eq. (14), we have, furthermore, 
employed the Markov approximation, which involves the shift of the integration limit J ' — > 
J °° and the replacement 



p(t-r)^e iCsT p(t). (16) 



The latter approximation is in line with the second order expansion in Hsb- Due to the 



shift in the integration limit, the master equation (14) is only valid for times longer than the 
correlation time of the bath |106f Tl08j. This is the case for the applications to steady-state 
transport considered in this paper, where only the long-time limit p(t — > oo) = p is required. 



Taken in the basis of eigenstates of the system Hamiltonian H$, the master equation (14) 
corresponds to the Redfield equation [9"9" | I1U3[ 109J. 

In the steady state transport regime, the above equation of motion becomes an algebraic 
set of equations 

= -i\H 3 ,p] (17) 
drtr B {[# S B, [H SB (T),pp B ]]}, 



OO 







which can be solved by standard linear algebra techniques. Thereby, the normalization 
constraint trg{p} = 1 ensures a unique solution. 



As discussed in the Introduction, the master equation (14) as well as its steady-state form 



Eq. (17), is valid for small molecule-lead coupling. Due to the neglect of terms of higher 
order in the system-reservoir coupling, it cannot describe tunneling in the non-resonant 
transport regime and misses the broadening of resonances due to this coupling. Except 
for these deficiencies, however, it provides a rather accurate description of vibrationally 
coupled electron transport in the resonant transport regime considered in this work. This has 
been demonstrated recently by comparison with results of nonequilibrium Green's function 
(NEGF) methods jH]. Test calculations show that all effects discussed in this work are 
also obtained with NEGF methods for vibrationally coupled resonant electron transport 

The basis functions, which we use to evaluate the reduced density matrix p and the 



master equation Eq. (17), are products of basis functions \a)\v) that span the subspace of the 
electronic \a) and the vibrational degrees of freedom \v), respectively. Thereby, the electronic 
basis functions are given in the occupation number representation i.e. \a) = \n1n2..), where 
rii G {0, 1} denotes the population of the ith electronic state. Throughout this article we 
consider a single vibrational mode with frequency Q. Hence, we represent the vibrational 
basis functions by harmonic oscillator basis functions where v G No stands for the 
excitation number of the vibrational mode. Thus, the coefficients of the reduced density 
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matrix can be written as 

ffi = {a\P*\J) = {amp\v>)\ul), (18) 

where upper case indices refer to states of the vibrational mode and lower case indices 
represent the electronic part of the respective Hilbert space. 



Evaluating Eq. (17) first between vibrational states (v\\ and |z/ 2 ), we obtain the equation 
-i{v x \ [H s ,p] h> = (19) 

kij V31/4 



with 

,1/1^2 

Here, denotes the Fermi distribution function of the respective lead, L or R, evaluated 
at energy e&, and S(x) stands for the Dirac-delta function, where e.g. 8(e + Uc\cx)\U) = 



5(e + 11). In Eq. (19), we have neglected all principal value terms. These terms describe 
the renormalization of the molecular energy levels due to the molecule-lead coupling |72j, 
which are negligible for the results discussed below. The thus obtained scheme represents a 
rate equation approach [7H IllOj . 



In the next step, Eq. (19) is evaluated with respect to the electronic basis functions. Since 



there is no further approximations involved, we present the rather lengthy expressions in 
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appendices [A] and |Bj where the results for a single electronic state and two electronic states 
can be found, respectively. 



C. Observables of interest 

1. Electronic population and vibrational excitation 

The diagonal elements of the density matrix, p™ a , encode the probability of finding the 
system in the product state |a)|z/). Hence, for a single electronic state on the molecular 
bridge the occupation of this state is given by the expression 

ni = {c\c x ) H = {c\c x )h (20) 
= tr s+B {^4 c i} = tr s {pc|ci} = J^PiJi- 

For two electronic states , the respective populations are given by 

ni = (c\c 1 )H = Y,PnM + Pil^ (21) 

V 
V 

Thereby, the subscript H/H denotes the Hamiltonian, which is used to evaluate the respec- 
tive expectation value. 

The average excitation of the vibrational mode involves a sum over all electronic degrees 
of freedom. For the transport scenario with a single electronic state on the molecular bridge, 
the average vibrational excitation reads 

A 2 

(a^a) H = (a^a) 7T + —(n 1 -25 1 n 1 + 5 1 ) (22) 
and respectively for the transport scenario with two electronic states 



(cta) H = V PZ + ^("i - + <*0 ( 23 ) 



A 2 

+ ^f(^2 - 2<5 2 n 2 + ^2 
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Since both operators c[, 2 Cy2 and era act in the subspace of the molecular bridge, the corre- 
sponding observables are fully determined by the reduced density matrix p. This is not the 
case for the current operator, which encompasses the bridge space and the subspace of the 
leads. 



2. Current 



The current through lead K } 1%, is determined by the number of electrons entering or 
leaving the lead per unit time (K G {L,R}) 



Ik = (Ik)h = -2e— ^(4c fc ) F 



(24) 



keK 



2ie 



VkMciXi)^ - *£<4xfo>2T 

_ ki ki 

Here, the constant (— e) denotes the electron charge and the factor 2 accounts for spin- 
degeneracy. The specific structure of the current operator requires the determination of 
Qq = (1 — P)g, since tis+B{P qIk} = 0. This projection of the full density matrix has 
already been used to derive Eq. ([13]) |104l 1105] , and reads 



Qg(t) 



-iQCt 



Qg(0) 



f dTe- iQCr QCPg{t-r) 
Jo 



(25) 



Using Eq. (25) to evaluate the expression for the current, Eq. (24), and employing the 



same approximations that were used for the derivation of the master equation, Eq. (17), the 
following expression for the current through lead K is obtained [62], [72} EH E2] 



IK 



—I 



drtr s+B { [#sb {t),pbp\ Ik}- 



(26) 



In the numerical calculations, Eq. (26) is further evaluated within the appropriate basis 
functions. The explicit formulas are given in appendices [A] and [Bj We note that the scheme 
described above is current conserving, i.e. 1^ = — Ir = /. 



D. Vibration in thermal equilibrium 



To analyze and identify vibrational nonequilibrium effects, it is often instructive to com- 



pare the results of Eqs. (17) and (26) with those, where the vibrational degree of freedom 
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of the molecular bridge is treated as a reservoir's degree of freedom. To this end, we define 
a new projection operator P' 



P'g(t) = pB'TMf?(*)} = PB'P'(t), (27) 

p B , = Z'-\- m *\ Z' = TrB>{e-^}, (28) 

where 

H B ' = ^e k c\c k + Sla)a. (29) 

k 

The thus defined reduced density matrix p' describes the electronic degrees of freedom on 



the molecular bridge. Analogously to Eq. (17) and following the same steps as outlined in 



Section II B we obtain the following master equation 



= -i[iW] (30) 
dr tr B / { [H S > B > , [H S > B > (r) , p p B >] ] } , 



with 

H 3 ' = ^li^Ci + ^Uij^Ci-SiX^Cj-Sj), 

i i<j 

(31) 

^s'B' = ^2(V ki X i clci + h.c.). (32) 

hi 

Treating the vibrational mode as a reservoir's degree of freedom does not exclude the pos- 
sibility to excite or deexcite it. However, it is assumed that such a nonequilibrium state 
relaxes on a short time scale to the thermal equilibrium state. This can be realized, e.g., by 
a strong coupling of the vibrational mode to a thermal bath. 

III. RESULTS 

We have applied the methodology outlined above to various models of vibrationally cou- 
pled electron transport in molecular junctions. The results described in this Section are 
structured according to the complexity of the models employed and the associated trans- 



port mechanisms. In Sec. HI A and B we present results for transport through a molecular 
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bridge with a single and two electronic states, respectively. Within each of these subsec- 
tions, we study first symmetric molecular junctions, proceeding with asymmetrically coupled 
junctions, and finally discuss, for transport through two electronic states, the influence of 
electron-electron interactions. 

In all cases considered the electronic states are coupled to a single vibrational mode 
with frequency Q = 0.1 eV. To represent the vibrational degree of freedom, the calculations 
employ iVb as = 200 vibrational basis functions, which provide converged results for all ob- 
servables and parameters considered in this section. As is shown in appendix [Cj vibrational 
coherences have no significant effect and are therefore neglected in these calculations. The 
temperature of the leads is set to 10 K in all calculations, since most experiments on molec- 
ular junctions are carried out at low temperatures. The same temperature is used for the 
vibrational mode in those calculations, where the vibration is treated in thermal equilibrium. 



A. Transport through a molecular junction with a single electronic state 

In this section we discuss vibronic effects in transport through a molecular junction that 
involves a single electronic state and a single vibrational mode. We start with a summary of 
the basic mechanisms of resonant emission and absorption processes (sketched in Figs, [l] and 
2]) for transport through a state that is symmetrically coupled to the leads. Furthermore, 
we discuss vibrationally induced rectification (23J 1111] and vibrationally induced negative 
differential resistance (NDR) gT\ US EH ESI EH M, M, E2] for junctions with asymmetric 
coupling to the leads. Understanding these generic mechanisms facilitates the discussion 



of transport through two electronic states in Sec. Ill B and extends our previous studies 
[4*0"] UU S3]- In particular, we show that the vibrational excitation of a molecular junction 
can only be understood if electron-hole pair creation processes are considered (cf. Fig. [2J [89- 
[91] . Furthermore, the results demonstrate that vibrationally induced rectification is strongly 
enhanced, if the vibrational degree of freedom, due to current-induced local heating, is in a 
highly excited nonequilibrium state. 
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1. Basic mechanisms 



We first consider a model with a single electronic state at the molecular bridge with 
energy ei = 0.6 eV and a moderate coupling to the vibrational mode, A = 0.06 eV. The left 
and the right leads are modelled as one-dimensional semi-infinite tight-binding chains with 
a semi-elliptic conduction band, for which the respective level- width functions T L / R (E) read 

Tl/r(E) = (33) 

l^l 2 \o, |£-// L/R |>2|/3|. 
Here, i^l/r = 0.1 eV denote the coupling strengths of the left and the right tight-binding chain 
to the electronic state at the molecular bridge and /3 = 3eV determines the band-width in 
both leads. The difference of the chemical potentials of the leads, — Hr = e$, determines 
the applied bias voltage $, which we assume to drop symmetrically at the contacts, i.e. 
= -/J-r = f$- 

Fig. [3] shows the corresponding current-voltage characteristics, the population of the 
electronic state at the molecular bridge (inset of Fig. |3]d) and the vibrational excitation of 
this system. To facilitate the discussion, the dashed line depicts results for a purely electronic 
calculation without vibronic coupling (A = 0). The corresponding current exhibits a single 
step at the bias voltage e$ = 2e±, indicating that for e$ > 2e\ electrons from the left lead 
can resonantly tunnel onto the electronic state and further to the right lead. These resonant 
transmission processes result in a current of I ~ 1.6 fih. 

The current- volt age characteristics and the vibrational excitation for a finite vibronic 
coupling strength of A = 0.06 eV are depicted by the solid black lines. For this case, the 
current rises in a multitude of steps. The first step appears at a lower bias voltage e$ = 2ej, 
reflecting the polaron shift of the electronic state. The following steps appear at voltages 
e$ = 2ei + 2nQ with n G N and a step height that gradually decreases with increasing 
bias voltage. The first step at e$ = 2l\ also marks the onset of resonant transport, where 
electrons from the left lead can resonantly tunnel onto the electronic state at the molec- 
ular bridge, accompanied by a transition from the initial vibrational state of the neutral 
molecule to the same vibrational state of the charged molecule. In a successive tunneling 
process from the molecular bridge to the right lead, which completes the transport process, 
the electronic-vibrational coupling may result in excitation (Fig. [l^,) or deexcitation (Fig. 
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FIG. 1: Basic schemes of vibrationally coupled electron transport processes for a single electronic 
state. Panels a) and c) depict examples for emission processes, where an electron sequentially 
tunnels from the left lead onto the molecule and further to the right lead, thereby singly exciting the 
vibrational mode of the molecular bridge (red wiggly line). Such emission processes are effectively 
'heating' the junction (local heating). An example for a respective absorption process is shown 
in Panel b), where an electron tunnels from the left to the right lead by absorbing a quantum of 
vibrational energy (blue wiggly line). These processes result in local cooling of the junction. 





FIG. 2: Example processes for electron-hole pair creation in a molecular junction. Panel a) 
depicts an electron-hole pair creation process with respect to the left lead by absorption of a single 
vibrational quantum. Panel b) represents an electron-hole pair creation process with respect to 
the right lead by absorption of two vibrational quanta. The absorption of two vibrational quanta 
normally occurs with lower probability. 

[Tjo) of the vibrational mode. At this bias voltage, there are mod (e$, CI) excitation (emis- 
sion) processes available (corresponding to the excitation of at most mod (e$, CI) vibrational 
quanta, as e.g. in Fig. [Lp), and a number of deexcitation (absorption) processes (Fig. |l]b) . 
For larger bias voltages, e$ > 2(ei + CI), electrons can excite the vibration upon tunneling 
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FIG. 3: Current and vibrational excitation for a generic model system with a single electronic state 
at the molecular bridge that is symmetrically coupled to the left and the right lead, and moderately 
coupled to a single vibrational mode. The inset shows the corresponding population n\ of the 
electronic level. The dashed line refers to a calculation where the electronic-vibrational coupling 
is set to zero. The solid gray and black line are obtained for an electronic-vibrational coupling 
strength A = 0.06 eV. Thereby, the gray line is calculated employing the thermal equilibrium state 
of the vibrational mode (at 10K), and the black line is obtained with its full current-induced 
nonequilibrium state. The solid red line depicts the vibrational excitation for this model system 
with a reduced electronic-vibrational coupling A = 0.03 eV. 
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from the left lead onto the molecule (Fig. [Tp), resulting in an increase of current and vi- 
brational excitation. The step-wise increase of the current associated with these processes 
gradually becomes smaller. This can be qualitatively rationalized by the Franck-Condon 
(FC) factors |Ao n | 2 = -7 (^"e"'^' 2 [73] that are associated with the transition proba- 
bility from the vibrational ground state to the nth excited state. For X/Q < 1, these FC 
factors decrease with increasing n. However, a quantitative description of the step-heights is 
more involved, since a variety of emission and absorption processes contribute to each step 
in the current- voltage characteristics [19J. 

We next consider the current-induced vibrational excitation depicted in Fig. [3]d. Similar 
to the current, the vibrational excitation increases in a step-wise way with increasing bias 
voltage. However, the step heights in the vibrational excitation become larger with increas- 
ing bias voltage. This is in striking contrast to the behavior of the current-voltage charac- 
teristics discussed above. Analogous to the current, the steps in the vibrational excitation 
are associated with the onset of resonant emission processes that involve successively more 
vibrational quanta. Hence, the relative step heights in vibrational excitation are expected 
to be larger than the relative step heights of the respective current-voltage characteristics. 
However, if only vibrational excitation/deexcitation processes induced by electron transport 
processes are taken into account, the resulting vibrational excitation should saturate as does 
the current. Moreover, for a smaller electronic- vibrational coupling a decrease of the rela- 
tive step heights would be expected due to the reduced Franck-Condon overlap of processes 
that involve multiple vibrational quanta. On the contrary, the comparison of the results for 
different vibronic coupling strengths in Fig. [3] shows the opposite behavior. 

A detailed analysis reveals that these intriguing findings are due to vibrationally in- 
duced electron-hole pair creation processes, schematically depicted in Fig. [2j The process 
of vibrational relaxation due to electron-hole pair creation is well known from spectroscopic 
[221 I113H115] and theoretical studies [551 - 1511 1116] of adsorbates at metal surfaces but has not 
been analyzed in detail in the context of nonequilibrium transport in molecular junctions 
[25| H2J S3]. In molecular junctions, such processes can deexcite the vibrational mode by 
creation of an electron-hole pair in one of the leads and contribute within the same order 
in the molecule-lead coupling V^i to the vibrational excitation as transport induced excita- 
tion/deexcitation processes {e.g. the one of Fig. [Tjo). If the bias voltage is increased, these 
processes are blocked one by one, because the resonant creation of an electron-hole pair 
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requires the absorption of increasingly more vibrational quanta. For typical values of A and 
the vibrational excitation, this means that the most important electron-hole pair creation 
processes are blocked first. As a result, the faster increase of the vibrational excitation with 
bias voltage is due to less efficient cooling by electron-hole pair creation processes |117j . 
Since the blocked electron-hole pair creation processes are even more important for smaller 
electronic-vibrational coupling, the associated step in the vibrational excitation characteris- 
tics becomes larger for smaller electronic-vibrational coupling. As a consequence, for large 
bias voltages the vibrational mode is more strongly excited the weaker it is coupled to the 
electronic state. This counter-intuitive behavior is related to the phenomenon of vibrational 
instability jJU IllOj . Hence, the observed increase of the step heights in vibrational exci- 
tation for increasing bias voltages as well as decreasing electronic-vibrational coupling is a 
result of the successive suppression of electron-hole pair creation processes. 

The current-induced vibrational excitation, in turn, has a pronounced effect on the 
current-voltage characteristics. This is demonstrated by the comparison of the solid black 
line in Fig. [3k with the solid gray line, which is obtained from a calculation with the vibra- 



tion kept in thermal equilibrium (cf. Section II D). In particular, the result for a thermally 
equilibrated vibrational degree of freedom shows a significantly larger current at a given 
bias voltage. For 2ei < e$ < 2(ei + Q), where (a^a) ~ 0.5, this behavior can be quali- 
tatively understood by considering only the vibrational ground (\v = 0)) and first excited 
state (\u = 1)). For the solid gray line, due to the low temperature T = 10K, the vi- 
brational mode is restricted to its ground-state, and the respective transition probability 
can be characterized by the FC factor |A"oo| 2 , which describes elastic tunneling processes 
from the left lead onto the molecular bridge. Note that at this bias voltage all relevant 
channels for tunneling from the molecule to the right lead are open, since ei — ep > 5fl. 
For the solid black line, the vibrational mode can be found in both states, and the re- 
spective transition probability can be characterized by a linear combination of FC factors: 
(l-(a t a))|A 00 | 2 + (a t a) (|An| 2 + |X i| 2 ). Thereby, the terms ~ |A 00 | 2 and ~ |An| 2 describe 
elastic tunneling processes from the left lead onto the molecular bridge, which is either in 
its vibrational ground- or first excited state. The third term ~ |Aoi| 2 represents inelastic 
tunneling processes, upon which the vibrational mode is deexcited. If X/Q < 1, such a linear 
combination of FC factors is smaller than |A o| 2 - Hence, for the given model parameters, 
vibrational excitation suppresses the first step in the current. Our results show that this 
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suppression of the current is a rather characteristic phenomenon, especially for larger bias 
voltages. This can be qualitatively understood with a similar analysis. If e.g. the bias voltage 
allows for m = mod (* — resonant emission processes with respect to the left lead (cf. 

Fig. ^p) and if the molecular bridge is in its ground-state (e.g. at T = 10 K), the transition 
probability from the left lead onto the bridge can be given as Yl™=o l^0n| 2 , which converges 
to unity with increasing m or bias voltage $. If the vibrational mode, however, is in a 
nonequilibrium state, where the population of the Zth vibrational level is given by cti ^ 8qi 
(I G N ), the corresponding transition probability is determined by J2i=o oo a i SLt^o l^«| 2 - 
Because vibrational excitation typically increases much faster than m does with increasing 
bias voltage, and because Y^n=o l^in] 2 ~ 1/2 for / ^> m, vibrational excitation typically 
results in a lower current. 



2. Vibrational rectification in asymmetric junctions 

In many experimental setups, single- molecule junctions are non-symmetric with respect to 
the left-right symmetry. This is the case in STM experiments [21H26] but often also in break- 
junction experiments [T2l [TH 152] . In this and the following section we study consequences of 
an asymmetric molecule-lead coupling on vibrationally coupled electron transport through 



a single electronic state. To this end, we employ the same model system as in Sec. Ill A 1 
but change the coupling of the electronic state to the right lead to z/r = 0.01 eV. 

Fig. [4] shows the respective current-voltage characteristics for different electronic- 
vibrational couplings A. If the vibration is kept in thermal equilibrium (Fig. Ek), the cor- 
responding current is approximately anti-symmetric with respect to bias, ~ —/(—$). 
Significant deviations from this antisymmetry appear around the onset of the current at 
e$ rs — 2e\. This can be understood in terms of tunneling processes at the boundary be- 
tween the molecule and the right lead, which represent the bottleneck for transport in this 
asymmetric transport scenario with i^r/^l = 0.1. Recall that the probability for tunneling 
processes, where an electron with energy e% enters the bridge from the right lead, is ~ |AToo| 2 , 
because the vibrational mode is essentially in its ground-state at T = 10 K. An electron with 
energy e\ + Q entering the bridge from the right lead will cause a single excitation of the 
vibrational degree of freedom. Because it is treated as a reservoir degree of freedom, the 
vibration will exhibit fast relaxation such that the next electron traversing the junction 
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finds the vibrational mode again in its ground-state. The respective transition probability 
is ~ \X i | 2 . Analogously, for electrons with energy Ii + nQ it is ~ |X 0n | 2 . Thus, for negative 
bias voltages, where electrons flow preferentially from right to left, the current increases with 
relative step heights that are determined by the FC factors |Xo„| 2 . In contrast to the findings 
in Sec. HI A 1 this is a quantitative statement. Thereby, the electronic state remains essen- 
tially unoccupied due to the asymmetry in the molecule-lead coupling, (c*c) ~ 2 f 2 ~ (cf. 
the inset of Fig. |4^i). For positive biases, on the other hand, the molecular electronic state 
becomes almost fully occupied, (ac) ~ v i\ v i ~ 1, as soon as the electronic state enters the 

R L 

bias window. Successive resonant emission processes with respect to the left lead are thus 
Pauli-blocked and do not result in further steps in the current- volt age characteristics. The 
step at e$ = 2ei, however, is already as high as the one in the electronic current without 
vibronic coupling and thus displays no suppression due to electronic-vibrational coupling. 
This is due to the fact that a rather large number h of resonant emission processes with 
respect to the right lead is already available at this bias voltage such that the sum of the 
respective FC factors approximately equals unity: ^2 i=0 fi \X 0i \ 2 « I. This requires the elec- 
tronic level to be located well above the Fermi energy such that e± — ep > nQ. Qualitative 
similar effects have been found previously in theoretical and experimental studies [231 HH] ■ 

Including the current-induced excitation of the vibrational mode in full nonequilibrium 
gives qualitatively different results (cf. Fig. [Ijo). In particular, for large negative bias voltages 
the current does not approach the maximum value given by the electronic current (A = 0), 
but remains at a significantly lower absolute value. Thus, vibronic coupling results in a 
persistent rectification of the current, which appears not only in the vicinity of e$ ~ — 2e± 
but for larger absolute values of the bias voltage as well. The comparison of Fig. |4^i and 
Fig. |4]3 shows that this vibrational rectification is a pure nonequilibrium effect. 

The respective vibrational excitation, depicted in Fig. [5] shows a similar asymmetry as 
the current. For negative bias voltage the level of vibrational excitation is much higher than 
for positive bias. We attribute this behavior to electron-hole pair creation processes with 
respect to the left lead (Fig. [2]). Due to the strong coupling of the left lead to the electronic 
state at the molecular bridge, these processes provide the most important cooling mechanism 
in the junction. For negative bias voltages, however, the creation of an electron-hole pair 
in the left lead requires the absorption of many vibrational quanta (> mod(2ei,f2) = 12 in 
the resonant transport regime), which is rather unlikely. For positive bias voltages, on the 



21 



< 



30 
20 
10 


-10 
-20 
-30 



< 
a 



: LL 



0.9 
0.6 

0.3 





-1 



n 



-l=0.12eV 
l=0.06eV 

-A=0.03eV 
- -A.=0eV 



-2-10 1 

bias voltage $ [V] 



30 
20 
10 

-10 
-20 
-30 



1 



0.9 
0.6 
0.3 





-1 




-X=0.12eV 
A=0.06eV 

-X=0.03eV 
- -X=0eV 



-2 -1 

bias voltage $ [V] 




FIG. 4: Current-voltage characteristics for a generic model system with a single electronic state 
that is asymmetrically coupled to the leads, and moderately coupled to a single vibrational mode 
with different coupling strengths A. The current-voltage characteristics of the upper panel have 
been obtained with the vibrational mode kept in thermal equilibrium, while for the ones in the 
lower panel the full nonequilibrium state of the vibrational mode is taken into account. The insets 
represent the corresponding populations of the electronic level. 

contrary, such an electron-hole pair may be generated by absorbing much less vibrational 
quanta, which is much more probable and effectively cooling the vibrational mode. This 
results in a lower level of vibrational excitation for positive bias voltages, but also in a 
much higher level of vibrational excitation for negative bias voltages. The enhancement 
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FIG. 5: Average vibrational excitation corresponding to the currents shown in Fig. HV 



of vibrational rectification, which we observe in Fig. |4j is a result of this higher level of 



vibrational excitation (cf. the discussion at the end of Sec. IIIAl). 



The phenomenon of vibrational rectification can be crucial for the interpretation of ex- 
perimental data. In particular, it explains the disappearance or suppression of vibrational 
side-peaks for a specific direction of the bias voltage, although such peaks may be clearly 
visible for the opposite direction of the applied bias voltage. This phenomenon has been 
observed in a number of experiments 



3. Vibrationally induced Negative Differential Resistance (NDR) 

In the model discussed in the previous section, the electronic state was located far away 
from the Fermi level of the system: 1\ — ep > 5f2. Thus the first step in the respective 
current- volt age characteristics involves a large number of emission and absorption processes 
that cannot be resolved separately. If the electronic state of such a molecular junction is 
located closer to the Fermi level, i.e. |ei — ep | < f2/2, every resonant emission (or absorption) 
process can be associated with a corresponding step in the transport characteristics. This 
is illustrated in Fig. [7| which represents the transport characteristics of a model system, 
where the energy of the molecular electronic state is e\ = 0.066 eV (ex = 0.03 eV). All other 



parameters are the same as in Sec. HI A 2 
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FIG. 6: Example transport processes that become active at a bias 
voltage 2(0 — ei), as the electronic state of our model system is 
located close to the Fermi-level of the system, i.e. \e\ — If\ < 0/2. 
If ei — ep > 0/2, the hopping process from the left lead onto 
the molecular bridge requires a higher bias voltage, <1? = 2ei > 
2(0 — ei), and in that case, these transport processes become active 
at the same bias voltage as e.g. elastic transport processes that do 
not change the vibrational state. 



We first consider results obtained with keeping the vibrational mode in thermal equilib- 
rium (at T = 10 K, i.e. essentially in its ground state), depicted by the blue and green lines 
in Fig. [7j Thereby, the blue line shows the current for positive bias voltages, /($), while the 
green line represents the current for negative bias voltages, —/(—$). Due to the asymmet- 
ric molecule-lead coupling, z/r/z/l = 0.1, distinct resonances appear whenever a tunneling 
process with respect to the right lead becomes available. This results, e.g. for the blue line, 
in a resonance at e$ = 2(Q — ei). The corresponding transport processes are schematically 
shown in Fig. [6j Similarly, resonances appear at e$ = 2\n{l — ei|, with n £ No, and relative 
step heights in the current that are associated with the FC factors |Xo n | 2 (as already out- 



lined in Sec. Ill A 2). Tunneling processes with respect to the right lead appear as steps in 



the green line as well, but at different bias voltages e$ = — 2(ei + nQ). As a consequence 
of the different position of the resonances, the number of tunneling processes with respect 
to the right lead differs for both lines at specific absolute values of the bias voltage, and 
therefore, the blue and the green line encircle small rectangles with a width of 4ei. The 
heights of these rectangles are given by the respective FC factors: |X 0n | 2 {n > 1). Outside 
these rectangles, the blue and the green lines have the same absolute values, corresponding 
to an equal number of active transport channels with respect to the right lead. 

If the vibrational mode is allowed to develop its full current-induced nonequilibrium state, 
we obtain the solid black line for positive, and the solid gray line for negative biases. Again, 
the black and the gray line encircle rectangles. However, the heights of these are enlarged 
compared to the ones encircled by the blue and the green lines. This is related to the fact 
that the current represented by the gray line decreases at bias voltages 2(ei — nfl), and 
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FIG. 7: Current and vibrational excitation for a generic model system with a single electronic 
state close to the Fermi-energy, |ei — ep| < 0/2, and a single vibrational mode. The dashed lines 
refer to a calculation with the vibrational mode in thermal equilibrium, while for the solid lines the 
vibrational mode is treated in its full current-induced nonequilibrium state. The cur rent- volt age 
characteristics for positive (black and blue line) and negative bias voltages (gray and green line) 
are overlayed to highlight areas, where /($) = — /(— <3?) holds. 

is thus significantly smaller than the one represented by the green line, while the black 
and the blue line display essentially the same current values at these bias voltages. For 
e<3> = 2(nQ + €i) (n > 1), when the black and the gray line start to overlap again, we observe 
a significant drop of the current represented by the black line. Interestingly, for both lines 
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such negative differential resistance coincides with the steps in the associated vibrational 
excitation characteristics, which is shown in Fig. [7Jd. Here, the black line represents the 
vibrational excitation induced by the current for positive bias voltages, while the gray line 
depicts the one for negative bias voltages. For both polarities of the bias voltage, the 
vibrational excitation increases at exactly those values of $, where electron-hole pair creation 
processes with respect to the left lead become blocked, i.e. for the black line at e$ = 
2{nQ + si) and for the gray line at e$ = 2(e x — nQ), where n > 1. Recall that electron-hole 
pair creation with respect to the left lead is the most important cooling mechanism due to 
the asymmetry in the molecule-lead coupling, = 10z/r. Weakening this cooling mechanism 
results in larger vibrational excitation, and consequently in a reduced current or negative 



differential resistance (cf. Sec. Ill A 1 ) . As a result, a vibrational nonequilibrium state may 



not only induce rectification, as pointed out in Sec. HI A 2 but also negative differential 
resistance (NDR). 

As is shown in Fig. [8j these features of negative differential resistance disappear one by 
one for stronger electronic- vibrational coupling A. The first one at e$ = 2(?i — Q), which is 
closest to the Fermi level, disappears for X/Q > 1 [731 US] . This can be understood, as before, 
by analysis of the tunneling processes with respect to the right lead. For bias voltages in the 
range — 2ei > e$ > 2(ei — Q) there is only a single electronic tunneling process available, 
which is associated with a transition probability |A o| 2 , because the vibration is not excited 
for this voltage. For smaller bias voltages in the range 2(ei — Q) > e$ > 2(ei — 2Q) resonant 
emission processes with respect to the left lead result in a finite vibrational excitation. As a 
consequence, the vibrational mode can be deexcited by a resonant absorption process with 
respect to the right lead. The corresponding transition probability for tunneling from the 
right lead onto the bridge thus involves a superposition of | A^oo | 2 and |An| 2 + |A"oi| 2 . For 
X/Q < 1, such a superposition gives a smaller transition probability than the previous one 
with only |A o| 2 , resulting in an overall smaller current (NDR). For X/Q > 1, this transition 
probability becomes larger than | A"oo | 2 • Consequently the current is also larger and the 
NDR at e$ = —2(Q — ei) disappears. The red line in the inset of Fig. [8] shows the difference 
|Aoo| 2 — (|An| 2 + |Aoi| 2 ) versus the electronic-vibrational coupling A. The transition from a 
smaller to a higher transition probability is represented by the zero-crossing at X/Q = 1. A 
similar analysis can be done for the second NDR feature at e$ = 2(li +Q), which disappears 
for X/Q > \/2. This behavior is depicted by the dashed green line in the inset of Fig. pj 
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FIG. 8: Current-voltage characteristics for a single electronic state, which is located close to the 
Fermi-energy, |ei — ep| < $7/2, and coupled to a single vibrational mode by different coupling 
strengths A. Features of NDR disappear one by one for increasing electronic-vibrational coupling. 
The inset shows the difference of the transition probabilities | ^Coo 1 2 — (|Xn| 2 + |Xoi| 2 ) ( re d line) 
and (|Xoo| 2 + |Xoi| 2 ) — (|Xn| 2 + |Xoi| 2 + |^12| 2 ) (dashed green line) as functions of the electronic- 
vibrational coupling strength A. The zero-crossings mark the value of the electronic-vibrational 
coupling, where the respective NDR- feature vanishes. 

which represents the difference (|X 00 | 2 + |X i| 2 ) - (|X n | 2 + |X 01 | 2 + |X 12 | 2 ). The NDR 
features at higher bias voltages vanish for increasing electronic-vibrational coupling as well. 
However, in this case a detailed analysis of all the contributing processes is more involved. 

Because this NDR effect is solely based on a vibrational nonequilibrium state, it is not 
exclusively restricted to asymmetric junctions. It is also found in symmetric junctions (see 
for example the gray line in Fig. 2 of Ref. HTjl . However, in symmetric junctions NDR 
is usually much less pronounced. This can be understood by the transport channels that 
become active, whenever electron-hole pair creation processes become blocked (e.g. the pair- 
creation processes in Fig. [2^, are blocked as soon as the transport channel represented by 
Fig. [TJ^ becomes active). In contrast to asymmetric junctions, all these transport channels 
contribute to the current and thus counteract the NDR mechanism that is induced by the 
higher level of vibrational excitation. 
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B. Transport through a molecular junction with two electronic states 



In this section we study transport through a molecular junction with two electronic states. 
In addition to the phenomena discussed above for transport through a single electronic state, 
three new aspects need to be considered jH]. First, in systems with multiple electronic states, 
the Hamiltonian in Eq. (|4]) comprises an electron-electron interaction term, which represents 
both Coulomb interactions and vibrationally induced electron-electron correlations. These 
electron-electron interactions result in a splitting of resonances jH] that depends on the 
specific population of the electronic states. For asymmetric junctions with a blocking state 
[15*1 [Ml [55] , they can induce strong NDR. While such NDR has been reported for specific 
values of the bias voltage $, here we propose another model system with a centrally localized 
electronic state, where NDR due to electronic-vibrational coupling extends over a broad 
range of bias voltages. Secondly, coherences of the density matrix can play a significant 
role in transport through multiple electronic states p"9j [80] . As is shown in appendix [Cj 
however, coherences are of importance only for asymmetrically coupled junctions with quasi- 
degenerate energy levels. For the systems considered in this Section, as in the preceding one, 
coherences of the density matrix can therefore be neglected. The third, and most intriguing 
aspect is that higher-lying electronic states facilitate the efficient absorption of vibrational 



energy [41, 42J. The role of these resonant absorption processes is elucidated in Sec. Ill B 1 



for symmetric molecular junctions, and in Sec. |III B 2 for asymmetric molecular junctions. 
This local cooling mechanism can be very efficient and crucial for the stability of the junction. 
Since closely-lying electronic states are typical for polyatomic molecules, this mechanism is 
expected to be of general importance in molecular junctions. It is also shown that repulsive 
Coulomb interactions may enhance this effect and thereby significantly improve the stability 
of a molecular junction ('Coulomb Cooling'). 



1. Resonant absorption processes via a higher-lying electronic state in symmetric molecular 
junctions 

We consider a model system with two electronic states, located at e\ = 0.15 eV and 
at €2 = 0.8 eV above the Fermi-level, respectively. The coupling strengths between the 
vibrational mode and the electronic states are given by Ai = 0.06 eV and A2 = — 0.06eV. 
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FIG. 9: Basic processes for vibrationally-coupled electron transport involving two electronic states, 
including resonant emission (a), resonant absorption (b) and electron- hole pair creation (c) pro- 
cesses. Due to resonant emission processes with respect to the lower-lying electronic state (panel 
a)), resonant absorption and electron- hole pair creation processes with respect to the second elec- 
tronic state become active even before this state enters the conduction window set by the applied 
bias voltage. 

In the results discussed in this section no Coulomb interactions are taken into account such 
that electron-electron interactions U12 = —2X1X2/Q are induced by vibronic coupling only. 
As before, the left and the right lead are represented by semi-elliptic conduction bands. The 
respective level-width functions T L / Rti j(E) thus read 



L/R,u 



(E) 



(34) 



^/r/l/r, j VlM^F, \E-»l/r\ < 21)81, 

\E-fi L/R \ >2\P\, 



0. 



where ^L/R,i = 0.1 eV denote the coupling strengths of the left and the right lead to the ith 
electronic state and /3 = 3 eV determines the band-width in both leads. 



Current-voltage characteristics for this model molecular junction are shown in Fig. [10 
Thereby, the dashed line represents results without electronic-vibrational coupling, A1/2 = 0. 
In this case, the two steps at e$ = 2ei and e$ = 2e2 are associated with the onset of 
transport through the lower-lying electronic state at €\ and the higher-lying electronic state 
at €2, respectively. 

The solid lines depict the results including the coupling of the two electronic states to the 
vibrational mode. We first consider the results (solid gray line), where the vibrational mode 
is kept in thermal equilibrium (T = 10 K). The respective current- volt age characteristics 
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FIG. 10: Current-voltage characteristics for a model molecular junction comprising two electronic 
states, which are coupled to a single vibrational mode. The dashed line corresponds to a calculation 
without electronic- vibrational coupling = 0. The solid black and gray line show results with 
a moderate coupling of the vibrational mode to both electronic states. For the solid gray line the 
vibrational mode is assumed to be in its thermal equilibrium state, while for the solid black line 
the full current-induced nonequilibrium state of the vibrational mode is taken into account. 

exhibits two major steps that correspond to the polaron shifted levels ei/2- Vibrational side- 
steps with respect to these resonances can be distinguished similar to the results obtained for 



a single electronic state in Sec. Ill A 1 (cf. Fig. [3J). A more detailed analysis shows, however, 
that the current obtained for the two-state system is not just the sum of the current through 
the individual states. For example, the step at e$ = 2e2 is just half as high as the one at 
e$ = 2?!. This behavior can be understood in terms of vibrationally induced electron- 
electron interaction, which is not present for a single electronic state. In this range of bias 
voltages, 2e2 < e$ < 2(t2 — 2\i\ 2 /fl), electrons in the left lead have not enough energy to 
doubly occupy the molecular bridge. Since the low-lying electronic state is half occupied 
at this bias, n\ = 1/2, the step at e$ = 2e 2 is thus reduced by a factor of 1 — n\ — 1/2. 
Consequently, this step can be associated with the electronically excited state of the anion. 
If the bias exceeds the value 2(e2 — 2\\\ij^t)^ electrons from the left lead can overcome 
the additional charging energy U . Hence, transport through the higher-lying electronic 
state becomes possible even though the low-lying electronic state is occupied, i.e. transport 
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through the di-anionic state of the junction. As a result, the single step at e$ = 2e 2 , which 
is associated with the onset of transport through the second electronic state, is split into 
two steps at e$ = 2e 2 and e$ = 2(e 2 — 2AiA 2 /f2). Vibrational side-steps with respect to 
the third, di-anionic resonance appear at bias voltages e$ = 2(e 2 — 2AiA 2 /f2 + nfl) (n EN). 
Note that U = — 2AiA 2 /f2 > since the vibrational coupling strengths Ai and A 2 differ by 
sign. For AiA 2 > the order of the steps associated with the excited state of the anion and 
the di-anionic state would be reversed. This scenario is described in Ref. HU 

Next, we consider the results obtained with the vibrational mode treated in nonequilib- 



rium (solid black lines of Figs. 10 and 11 ). The current-induced excitation of the vibrational 
mode changes the current- volt age characteristics profoundly. The most striking difference 
is the rise in current even before the second electronic state enters the bias window. This 
rise of the current is facilitated by the absorption of one or more vibrational quanta (see 
Fig. [9|d for an example process). These resonant absorption processes appear at voltages 
e$ = 2(e 2 — nfl) and e$ = 2(e 2 — 2AiA 2 /f2 — nfl), and can take place only if the vibrational 
mode is in an excited state. Resonant emission processes with respect to the lower-lying 
electronic state (cf. Fig. [9k), however, are efficiently exciting the vibrational mode, and 



thus, provide the vibrational energy required for these processes (see Fig. 11). In addition, 
these resonant absorption processes result in a pronounced broadening of the resonances 
that are associated with the second electronic state. This leads to an almost Ohmic con- 



ductance characteristics, which is observed, e.g. in Fig. 10 This broadening complicates the 
spectroscopy of molecular levels in single-molecule junctions \§T\ 192"] - 

Resonant asorption processes due to higher-lying electronic states may have an even 
more profound effect on the vibrational excitation of the molecular junction, as shown in 
Fig. [TTj For lower bias voltages, $ < 0.8 V, the vibrational excitation depicted by the black 



line (which was obtained for the same parameters as the curves in Fig. 10) exhibits 
an increase similar to the case of a single electronic state (cf. Fig. ^>). For larger bias 
voltages, however, the vibrational excitation drops by more than 50%, before it starts to 
increase again for e$ > 2e 2 . This pronounced reduction of vibrational excitation is caused 
by absorption of vibrational energy via resonant absorption processes with respect to the 
higher-lying electronic state. Since the second rise of vibrational excitation is shifted by 
more than IV, the molecular junction is effectively stabilized over a wide range of bias 
voltages [HI H2] . It is noted that a similar decrease of vibrational excitation with increasing 
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FIG. 11: Vibrational excitation for the two-state model molecular junction employed for the 



current- voltage characteristics of Fig. 10 (black line). The different results correspond to different 



energies €2 for the higher-lying electronic state. The solid black line shows the vibrational excitation 



that corresponds to the respective cur rent- volt age characteristics of Fig. 10 The solid blue and 
red line depict results, where the higher- lying electronic state is located closer to the lower- lying 
electronic level. 

bias voltage for a molecular junction has been observed in recent experiments by Ioffe et al. 

The details of this stabilization or cooling mechanism depend on the energy gap between 
the higher- and lower lying electronic state. In particular, if the higher-lying state is located 



too close to the lower- lying state (see the solid red and blue line in Fig. 11 ), resonant emission 
processes with respect to both states become active at the same time, and no decrease of 
vibrational excitation is observed. Nevertheless, we expect this cooling mechanism to be 
relevant for most molecular junctions, because polyatomic molecules often exhibit multiple 
closely-lying electronic states. 
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2. Resonant absorption processes via a higher-lying electronic state in asymmetric molecular 
junctions 



Resonant absorption processes with respect to a higher-lying electronic state involve 
not only electron transport processes (as in Fig. [9b) but also electron-hole pair creation 



processes (sketched in Fig. |9p). As was discussed in Sec. Ill A, electron-hole pair creation 
plays an important role especially for asymmetric junctions, where the electronic states on 
the molecular bridge are coupled to the leads with different coupling strengths Vk,%- There 
are eight topologically different scenarios for coupling two electronic states asymmetrically to 
a left and a right lead. Since all scenarios show similar effects, we focus in this section on the 



physically most relevant, where the two states of the model system introduced in Sec. Ill B I 



are strongly coupled to the left, but weakly coupled to the right lead, ^,1/2 — 10^r,i/2 = 
0.1 eV. This coupling scenario describes, e.g., the experimental setup of a scanning tunneling 
microscope (STM) [221I21H2S]. 

The corresponding current- volt age characteristics of this model is represented by the 



solid red line in Fig. [12p ,. For comparison, we also show the current of the corresponding 
symmetric model molecular junction with ^l,i/2 — u k,i/2 = 0.1 eV (solid black line) rescaled 
by a factor of 1/100. For positive bias voltages, the current shows only two steps at e$ = 2l\ 
and e$ = 2(t2 — 2X1X2/0). There is no splitting of the resonances associated with the higher- 
lying electronic state due to electron-electron interactions U, since the low-lying electronic 



state is almost fully occupied, ri\ = 1 (cf. the inset of Fig. 12 1), once it enters the bias 
window. Transport through the electronically excited state of the anion is thus not visible in 



the respective current- voltage characteristics (cf. the discussion in Sec. Ill B 1). For negative 
biases, a multitude of vibronic resonances are seen resulting in the same rectifying behavior 



that was already discussed for a single electronic state in Sec. |III A 2[ Small NDR features 



appear as well (highlighted by red arrows in Fig. 12), and can be related to the same 



mechanisms that was discussed in Sec. HI A3 For positive bias voltages resonant absorption 
via the higher- lying electronic state does not significantly contribute to the current, since the 
vibrational mode is efficiently cooled by electron-hole pair creation processes with respect 
to the left lead via both electronic states (see, e.g., Fig. [9p). These processes are dominant 
in this regime due to the asymmetry in the molecule-lead couplings vk,%- For negative 
biases, however, electron-hole pair creation with respect to the left lead becomes inefficient 
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FIG. 12: Current-voltage characteristics and vibrational excitation for a molecular junction with 
two electronic states that are moderately coupled to a single vibrational mode. The solid red line 
refers to a calculation with asymmetric molecule-lead coupling, where both electronic states are 
strongly coupled to the left lead and weakly coupled to the right lead: ^l,i/2 = 10^r,i/2- The solid 
black lines represent the result for the corresponding symmetric junction, which are the same as 



the solid black lines in Figs. 10 and 11 The corresponding current-voltage characteristic is thereby 
rescaled by a factor of 1/100 for a better comparison with the red line. The inset in the upper 
panel shows the population of the electronic states n%/2 for the asymmetric junction. 
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due to the increased amount of vibrational energy required for these processes. Therefore, 
vibrational excitation increases faster for negative bias voltages than for a symmetrically 
coupled molecular junction. This asymmetry in vibrational excitation is depicted in Fig. 
[l2p . Hence, in contrast to our previous findings for a symmetric junction, where the steps 
associated with the higher-lying electronic state are significantly broadened, we conclude 
that spectroscopy of molecular orbital energies may be more easily performed with a STM- 
like setup. Although resonant absorption processes are active, the mechanism of vibrational 
rectification restores the signatures of the individual molecular states. 



3. Coulomb cooling 

In this and the following subsection, we study the effects of repulsive Coulomb interactions 
(U > 0) on vibrationally coupled electron transport through a molecular junction. We first 
consider the influence of these interactions on resonant absorption processes, and thus on the 
stability of such a junction. NDR effects that arise due to repulsive Coulomb interactions 
will be considered in Sec. IHI B 41 



To this end, we employ the same two-state model system as in Sec. HI B 1 , but add an 
additional charging energy of U = 0.5 eV in the Hamiltonian, Eq. ([T]). This term accounts for 
repulsive Coulomb interactions between two electrons that are occupying the two electronic 
states of the junction. As a consequence, the resonance associated with the di-anionic state 
is shifted towards higher energies, from I2 — 2A A A2 to I2 — 2A A A2 + U . The resulting current- 



voltage characteristics (data not shown) shows the same features as analyzed in Sec. Ill B 1 
The major difference between the two scenarios is a shift of all steps associated with the di- 
anionic state towards higher bias voltages. As a result, the current increases more slowly in 
the resonant transport regime for e$ > 2e\. The respective vibrational excitation is shown in 



Fig. 13 1. Thereby, the solid black, red and blue lines directly correspond to the lines in Fig. 
11 which were obtained using the same model parameters but without Coulomb interaction. 
Since resonant absorption processes with respect to the di-anionic resonance are shifted to 
higher energies, the level of vibrational excitation increases faster than without repulsive 
Coulomb interactions for lower bias voltages. At higher bias voltages, however, these cooling 
processes are more efficient than without Coulomb interaction, as they require the absorption 
of less vibrational quanta. Moreover, they are not competing with absorption processes 
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FIG. 13: Upper Panel: Vibrational excitation of a model molecular junction with two electronic 
states, both of which are moderately coupled to a single vibrational mode and symmetrically to 
the leads. In addition, a repulsive electron-electron interaction U = 0.5 eV is taken into account. 



Lower Panel: Difference between the vibrational excitation shown in Fig. 11 and the upper panel, 



i.e. with and without Coulomb interaction, (a^a)H. 



{a)a) H , 



C/=0.5eV- 



associated with the excited state of the anion. As a result, the vibrational excitation, which 
we obtain for large bias voltages including a repulsive Coulomb interaction U = 0.5 eV, is 
reduced compared to that without additional electron-electron interactions. Thus, repulsive 
Coulomb interactions effectively stabilize a molecular junction. This cooling mechanism, 



here referred to as 'Coulomb Cooling' is analyzed in more detail in Fig. 13 d. The three 
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different lines depict the difference in vibrational excitation, (a^a)n,u=o ~ ( aJ *a)H,u=o.5eV, 
obtained without, U = OeV, and with Coulomb interaction, U = 0.5 eV, respectively. The 
comparison shows that, except for small bias voltages, the Coulomb interactions lead to a 
strong reduction of the vibrational excitation and thus significant cooling |118j . 



4- Negative Differential Resistance (NDR) induced by Coulomb interactions 

In this last section, we study NDR effects induced by Coulomb interaction. 

First, we consider junctions with asymmetric molecule-lead coupling. As we have already 
seen, asymmetries in the molecule-lead couplings can cause strong changes in the popula- 
tions of the electronic levels. This is particularly important in the presence of Coulomb 
interactions, where they may result in pronounced NDR effects. To analyze this effect, we 
employ the same model parameters as in the previous section, except that the coupling of 
the higher-lying electronic state to the right lead is reduced to z/a^ — 0.01 eV. This param- 
eter regime is known as a blocking-state scenario [4"5l EU |95]. Furthermore, the energy of 
the higher-lying electronic state is chosen as e 2 = 0.4 eV. The current- voltage characteris- 



tics depicted in Fig. 14 shows that the combination of an asymmetry in the molecule-lead 
coupling and Coulomb interactions results in a significant NDR effect at a bias voltage of 
e$ = 2e 2 . Once the bias exceeds this value, e<I> > 2e 2 , the second electronic state becomes 



almost fully occupied (see the inset of Fig. 14). This effectively blocks electron transport 



through the first electronic state, because an additional charging energy U is required, and 
the respective current drastically decreases. If the bias is increased further e<3? > 2(ei + U), 
transport through the first electronic state can take place, although the second electronic 
state remains partially occupied. For negative bias voltages, however, the second electronic 
state leaves almost no traces in the current- volt age characteristics, as it remains unoccupied. 
As a consequence, NDR is not observed for this direction of the bias voltage. 

Another interesting mechanism for NDR is observed if the higher-lying electronic state of 
the molecular junction is weakly, but symmetrically coupled to both leads, z/ L / r; 2 = 0.01 eV. 
Such a state corresponds to a molecular orbital, which is stronger localized in the central 



part of the junction. The solid red line in Fig. 15 represents the corresponding current- 
voltage characteristics. In contrast to the NDR effect discussed above, the current decreases 
over a broad range of bias voltages, from 2ei to 2(ei + U). Even for negative bias voltages, 



37 



< 

3L 



CD 



Q 



-1 




0.9 






0.6 






0.3 











t 








// 

NDR 



symmetric junction 

blocking state 



1 1 

bias voltage <J? [V] 



FIG. 14: Current- voltage characteristics for a two-state molecular junction with a single vibrational 
mode, where a repulsive electron-electron interaction of U = 0.5 eV is taken into account. The solid 
black line represents the current for a symmetrically coupled junction, and the solid red line the 
current for a junction with a blocking state |45|, IMj I95| . i.e. a higher- lying state that is weakly 
coupled to the right lead. The populations of the electronic states which correspond to the 
current- voltage characteristics depicted by the solid red line, are shown in the inset. Therein, the 
solid black line displays the population of the blocking state. 

the absolute value of the current decreases from — 2ei to —2(ei + U). Hence, this mechanism 
for NDR is symmetric with respect to bias. The decrease in the current is accompanied by 
an increase in the population of the second electronic state, which, as before, is blocking 



the current through the first electronic state. The inset of Fig. 15 depicts the populations 
ni/2 of both electronic states and shows the successive increase of n%. Thereby, the higher- 
lying state, due to resonant absorption processes, becomes occupied even before it enters 
the bias window. Increasing the bias voltage, the occupation number n 2 increases even 
further, as electronic tunneling processes and resonant emission processes become active. 
Overall, this increase takes place in small steps, resulting in an almost smooth current-voltage 
characteristics. Thus, transport through the lower-lying electronic state is not blocked in 
a single step, as for the NDR discussed before, but in a succession of several small steps, 
which results in a NDR feature that extends over a broad range of bias voltages. 
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FIG. 15: Current-voltage characteristics for a symmetrically coupled molecular junction, where 
two electronic states are moderately coupled to a single vibrational mode and exhibit repulsive 



Coulomb interaction, U = 0.5 eV. The solid black line, as in Fig. 14, represents the result, where 
both states are coupled to the leads with the same coupling strengths vk,%- For the solid red line 
the coupling of the higher-lying electronic state to both leads was decreased to ^l/r,2 = O-I^l/r.,1- 

IV. CONCLUSION 



In this paper we have studied vibrationally coupled electron transport in single-molecule 
junctions. The study was based on generic models for molecular junctions, which include 
electronic states on the molecular bridge that are vibrationally coupled and exhibit Coulomb 
interaction. The transport calculations employed a master equation approach, which is 
based on a second order expansion in the molecule-lead coupling. To the given order in the 
molecule-lead coupling, electronic-vibrational coupling and Hubbard-like electron-electron 
interactions are accurately described in this approach. 

The results obtained for a series of models with increasing complexity, show a multitude 
of interesting transport phenomena, including vibrational excitation, rectification, negative 
differential resistance (NDR) as well as local cooling. While some of these phenomena have 
been observed or proposed before, the present analysis extends previous studies and allows 
a more detailed understanding of the underlying transport mechanisms. In particular, the 
analysis shows that many of the observed phenomena cannot be explained if only transport- 
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induced processes (cf. Fig. [T]) are taken into account, but also require the consideration of 
electron-hole pair creation processes (cf. Fig. [2]). For example, these processes explain the 
increase of the step heights in vibrational excitation with increasing bias voltage and/or 
smaller electronic-vibrational coupling. Our results also show that in junctions with asym- 
metric molecule-lead coupling, electronic-vibrational interaction may result in a significant 
rectification of the current and the vibrational excitation, and even negative differential re- 
sistance. Both phenomena are accompanied by an highly excited nonequilibrium state of the 
vibrational mode. In these cases, electron-hole pair creation processes play a key-role. In 
contrast to transport-induced processes, they involve only one of the leads, and thus, effec- 
tively transfer the asymmetry of the molecule-lead coupling to the transport characteristics. 

Extending previous work |4T] , we have also given a detailed analysis of the influence of 
multiple electronic states on vibronic effects in molecular junctions. The results show that 
resonant absorption processes involving higher-lying electronic states efficiently reduce the 
level of vibrational excitation and thus can stabilize a molecular junction over a broad range 
of bias voltages jmH2]. In this context, we have shown that repulsive Coulomb interactions, 
which shift these resonant absorption processes to higher energies, can strongly enhance this 
cooling mechanism ('Coulomb Cooling') and thus improve the stability of the junction. 
Repulsive Coulomb interaction may also cause pronounced negative differential resistance 
(NDR). While such NDR effects were reported before at specific values of the bias voltage 
EH ES], in the model proposed in this work NDR is caused by vibronic coupling and 
extents over a broad range of bias voltages. It should be noted that single- molecule junctions 
often exhibit a number of closely-lying electronic states, which are strongly coupled with the 
junctions vibrational degrees of freedom. Thus, it can be expected that the phenomena 
analyzed in this work are of relevance for many molecular junctions and may facilitate the 
understanding and interpretation of experimental data. 
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Appendix A: Master equation and current formula for a single electronic state 



In this appendix, we give the explicit formulas for the master equation, Eq. (19), and 



the current, Eq. (26), for transport through a single electronic state coupled to a sin- 
gle vibrational mode. The density matrix of this scenario consists of 4iV^ as elements, 
PcmT j Pio^ 2 > tfoi 2 > Pi 1 ! 2 > corresponding to the two charge states |0) and In Eqs. (19) 



and (26) coherences between different charge states p^J 2 and p^ 2 are not coupled with the 
elements Pq 1 ^ 2 and p" 1 ^ 2 that are diagonal with respect to the electronic degrees of freedom. 
Therefore, these elements do not need to be considered in the following. If we evaluate Eq. 



(19) in between the charge states (0| and |0), we obtain 

-2m{vx - u 2 )p u ^ 2 = (Al) 
/] Ik (E U3U4 )T k (E U3U4 )X UlU3 Xl aUi pQ^ 2 



U3U4K 

- /j (1 ~ fK(E U4U2 ))r K (E U4U2 )X UlU3 Xl 4 „ 2 p l [ 3 J 4 



V3V4K 



/ y (1 - fK(E U3Ul ))T K (E U3Ul )X UlU3 Xl 4V2 p^J 4 



+ ^2 ft<{Eu A u 3 )^K{E UA y 3 )X U3UA Xl iU2 pQ\ 
V3V4K 



) 



with 



E Ua u b = ei + Q.{v a - u b ). 
The corresponding (1|..|1) component reads 

-2ift(>i - v 2 )p\^ 2 = (A2) 
^ (1 - fK{E V4V3 ))T K {E V4V3 )X V3V4 Xl xVz p v ^ 2 

- 22 f K ( E ^4)^K(E u . 2U4 )X U4U2 Xl iUa p^Q 4 

u 3 i/ 4 K 

- 2J fK(E UlV3 )T K (E UlU3 )X U4U2 Xl iU3 p^Q 4 
U3U4K 

+ ^ (1 - fK(E U3U4 ))T K^E^^X^Xl^p^J 3 . 
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From this set of algebraic equations, Eqs. (Al) and (A2), we determine the elements of 



the reduced density matrix, which are used to compute the current-voltage characteristics 
and the corresponding vibrational excitation of a single-molecule junction represented by a 



single electronic state (cf. Sec. II C). Thereby, the current through a single electronic state 
is explicitly given by 



IK 



(A3) 



+ e fK{E UlU2 )T K (E UlU2 )X UZUl Xl iV2 p u 2 Q 

+e ^2 fK(Eu 3 U2)^K(Eu 3 u 2 )Xu2U3 X l 3 u 1 Po 1 ,0 

-e 2j ( l - Sk{E V2V3 )) T K (E U2U3 )X U3Ul Xl 2U3 p^l 



>>':>, 



V1U2V3 



In Eq. (A3), principal value terms, as for the computation of the reduced density matrix, 



are disregarded. 



Appendix B: Master equation and current formula for two electronic states 



In this appendix, we give the explicit formulas for the master equation, Eq. (19), and the 



current, Eq. (26), for transport through two electronic states coupled to a single vibrational 



mode. Similar to the scenario with one electronic state, coherences between different charge 
states are not considered, as they are not coupled with the remaining elements of the density 
matrix. Thus, we are left with the computation of six different density matrix elements: 
Po7oo> PT (hJ e {1, 2}), and P 1\% where we employ the states |00), |1> = |10), |2> = |01), 
and 1 11). For notational convenience we introduce indices i, defined by 1 = 2 and 2 = 1. 
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Accordingly, we evaluate Eq. (19) between the states (00 1 and 1 00) 



anfa - z/ 2 )p Xo 2 o = (bi) 



iju s i>4,K 
ij 1/31/4 K 

+ ^ K > ii (Eiv 4 u 3 )fK(E iu4U3 )Xi jU3U4 X} U4U2 P QQ , 
W3V4K 



with 



E lVaUh = €i-U 126 ■ + tt(v a - v h ). 



Evaluation of Eq. (19) between the states (o| and \b), where a, b G {|10), |01)}, gives 

- 2t (n(u x - v 2 ) + e a - t b ) p v a f = (B2) 

^ ^ K,ia(E^ u . iU4 ) fl<(E+ U3U J Xa^ vz X\ V3V4 P^ 2 (-l) a+l 

+ XI ^aj^^^X 1 _ fK{Ej v ^)x i ^ VA x'i vxv ^p 

^3^4 / 1 \a+b 



"~ ^ ^ K > ab ^b,V2V4> 'J F K(E bu2U JXb tUi v 2 Xl UlU3 P QQ ( 
V3V4K 

- ^2 ^ Kfia{ E i^ iUl ){l - fK(E^ U3Ui ))Xa, UlU3 X^ 
1/3U4 K 

•Plfl 4 l(-l) a 



^ K,ab(E aUiU3 ) fK(E aViU3 )Xb )U4U2 Xl UiU3 P QQ Q 

V3V4K 

+ X ^^(^S^X 1 ~~ fKiE^vJJX^^Xl^piy 
IV3V4K 



with 



= ei + ^l-^+QK-!/,). 
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The respective (11|..|11) component reads 



- Unfa - u 2 )p^ = (B3) 
^ Fk,u(E£ U4U3 )(1 - fK(E^ U4U .J)X i:U . iV4 Xj jViU3 p^ 1 

iv 3 v^K 

~ E rK,ij(E£^ v Jf K (E£^ v jX j>V MXl^pg v *(-i) t+3 

ijv 3 v 4 K 



- E T K,ij(E^ UiV3 )f K (E^ UiU jX jtU4U2 Xl UiU3 p^ U4 (-iy +: > 

+ E r ^(^3^)( 1 _ fK{E^v z vJ) X i,VAV2 X \,v z V4P V ll,lV 



iv 3 v^K 

As for a single electronic state, this set of algebraic equations determines the density matrix 
describing transport through two electronic states. The current through two electronic states 
is explicitly given by 

1 



-Ik = (B4) 

e 

- ^2 ( 1 ~~ fK(E itV2Ul ))T K ji(E ijV2Ul )Xi tUlU2 Xl U3Ui p^ 3 



~ ^2 ( l ~ f K{Et >V2Ul )) ^K,ii{Ef tU2Ul )X i ^ lU2 X\ V3Ui p u ^l l 

i,v\,v 2 ,v 3 

+ E fK(E^ VlU2 )T Ktij (E^ UiU jX jtU3Ul Xl UiV2 p^ 

+ ^2 fK(Ei tUl u 2 )^K,ii(Ei tUlU2 )Xi :U3Ul Xj UiU2 P l 
i,u 3 u 2 )^ K,ii(Ei tUa u 2 )Xi^ U2U3 X iu3Ui PQQ 



2^3 
1,00 



1^2 

1,00 



,V 2 ,V 3 



_ E ( X ~ f K ( E i,V2V 3 )) ^K,ii(E^u 2 u 3 ) X i,u 3 u 1 Xj tV2U3 p'^ 2 
i,v\,v 2 ,v 3 



y-i.) JK yi^ i v3U2 )L K jiyi^ i u3U2 )SY i)P2V3 X ] 
1 _ f.,(TP+ \\r.,..(~p+ W. V 

:,iaj ^3^11,11 

- E ( 1 ~ ^ (Ei,u 2 u 3 ) ) rv^j (E ijV2U3 )X j>V3Ul Xj^^pf^ 2 . 

Thereby, principal value terms, as for the computation of the reduced density matrix, are 
disregarded. 
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Appendix C: Effects of electronic and vibrational coherences 



None of the results that we have discussed in Sec. |III| are significantly influenced by 
electronic or vibrational coherences of the reduced density matrix. This is due to the specific 
model parameters used. The respective eigenstates \a)\u) (a e {0, 1} or a G {00, 01, 10, 11}) 
do not exhibit any (quasi-)degeneracies. However, in realistic systems the molecular states 
\a)\v) can show some (quasi-) degeneracy, and therefore, we investigate the effect of such 
coherences in this section. 

For a single electronic state and a single harmonic mode, the molecular states are quasi- 
degenerate, if the broadening of the levels exceeds the level spacing set by the frequency fl. 
But in this regime, where Q < T, a perturbative expansion in V& may not be appropriate. 
Therefore, we start our discussion with a model system comprising two electronic states. 
We expect coherences to play a major role for degenerate or quasi-degenerate electronic 
states, i.e. for \ei — e 2 | < T [72\. As long as the couplings of the two quasi-degenerate 
electronic states to the leads are symmetric, we find the same results with and without 



coherences. This can be understood by inspection of Eqs. (Bl) and (B3). For symmetric 



junctions with €\ ~ e 2 and with U ~ 0, where poo,oo ~ Pii.n holds, Eqs. (Bl) and (B3) 
have almost the same structure. The only difference, in that case, is the sign by which 
coherences enter these equations. For that reason, the coherences pi j2 = p 2 i cancel in these 
equations and, as a result, do not influence the respective transport characteristics. Only 
for non-symmetric couplings to the leads we find a significant effect of coherences on the 
transport characteristics of this system. The most pronounced effect appears, if one of the 
molecule-lead couplings differs by sign e.g. for ^l,i/2 = z/r i = — z/r j2 and if the two states 
are degenerate e 2 = £\. This specific model system can be identically transformed to two 
orthogonal states, which are not interacting with each other and which are coupled to one 
of the leads only (either left or right). Hence, the current through this system is zero for 
any bias voltage W7\. However, if we disregard electronic coherences, we obtain a finite 
current that corresponds to the current of two states symmetrically coupled to the leads. 
Thus, for quasi-degenerate electronic states that are non-symmetrically coupled to the leads, 
electronic coherences must be accounted for to obtain physically correct results. 

The role of vibrational coherences can be studied employing a similar model system with 
two electronic states and energies that differ by the frequency of the vibrational mode, i.e. 
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^2 = ei+fi. Systems with €2 = ti+nVt and n > 2 display similar but less pronounced effects, 
as the impact of coherences decreases the further away they are located from the diagonal 
of the density matrix [72]. Again, for totally symmetric molecule-lead couplings vk,%-, we do 
not observe a significant influence of coherences. Only for asymmetric transport scenarios 
we find coherences to play a significant role for the transport characteristics. In Fig. 16 we 



present the current- voltage characteristics and the vibrational excitation for the asymmetric 



model system introduced in Sec. |III B2] with the energy of the higher-lying state adjusted to 
€-2 = t\+Vt = 0.25 eV. Thereby, the red line represents a calculation where all coherences are 
disregarded, while the green line represents the results for a calculation where all coherences 
are taken into account. Only in the vicinity e$ = 2e\ to e$ = 2(e~2+U) vibrational coherences 
influence the current-voltage characteristics. For positive bias voltages, the first step in the 
green line is diminished, since coherences result in a small population of the second higher- 
lying electronic state (cf. the inset of Fig. [T6^l) , which is thus blocking transport through the 
low-lying electronic state due to vibrationally induced repulsive electron-electron interactions 
U = —2X1X2/^1. This blocking is lifted again for higher bias voltages, when electrons in the 
left lead have enough energy, i. e. more than Ii+U. Vibrational excitation is enhanced as well, 
because there is an additional resonant emission process for tunneling from the higher-lying 
electronic state to the right lead. Similarly, coherences result in a somewhat larger current 
for bias voltages 2e2 < e$ < 2(6% + U), where the higher- lying electronic state enters the 
bias window. Again, the population of the second electronic state is increased by vibrational 
coherences, although repulsive electron-electron interactions U block the population of, and 
similarly, transport through this state. Coherences soften this blocking, resulting in a larger 
current and vibrational excitation. 

We expect systems with more than one vibrational degree of freedom to display gradually 
more quasi-degenerate levels, and hence, coherences to play a gradually more important 
role. Moreover for anharmonic potentials that are describing e.g. molecular motors |119j . 
coherences are crucial to characterize the actual motion of the molecule. 
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levels, where the dashed lines refer to the population of state 1 (m) and the solid lines to the one 
of state 2 (712). 
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